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Determining local transientness of audio signals 

Stephane Molla and Bruno Torresani 



Abstract 

We describe a new method for estimating the degree of "transientness" and "tonality" of a class of compound 
signals involving simultaneously transient and harmonic features. The key assumption is that both transient and 
tonal layers admit sparse expansions, respectively in wavelet and local cosine bases. The estimation is performed 
using particular form of entropy (or theoretical dimension) functions. We provide theoretical estimates on the 
behavior of the proposed estimators, as well as numerical simulations. Audio signal coding provides a natural field 
of application. 

Index Terms 

audiophonic signal, transient, tonal, wavelet basis, local cosine basis, sparsity. 
EDICS: l.TFSR, 2.AUEA 

I. Introduction 

Many generic signal classes feature significantly different "components", such as transients, (locally) 

sinusoidal or harmonic "partials", or stochastic-like components in sounds, or edges, textures, etc. in 

> 

k> ' i ma g es - Detecting the presence of such components is one of the classical signal processing problems 

(see for example [1] and [2] and references therein for reviews.) Another interesting problem is to estimate 
whether a given portion of a signal is for example more transient than harmonic or periodic, or in other 
words to estimate "transientness" or "tonality" indices: quantitative measures of the local proportion of 
transient and tonal features in a signal. Such indices find immediate applications in several contexts, 
including the hybrid signal coders [O, [@], ED which use different methods for encoding transient or 
tonal regions (and were the main motivation of this work), more general purpose hybrid models Q5], or 
similar recent ideas in image coding [J6l, flV). While there exist fairly standard tools for transient detection 
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or harmonic signal detection, the problem of quantitative measure of proportion does not seem to have 
received much attention. 

We propose here simple criteria, based on transform coding ideas, for estimating such indices. The main 
idea is to use orthonormal bases in signal spaces which are significantly different from each other in the 
following sense: a given component has a sparse expansion in a given basis, while the others have dense 
expansions. Information theoretic criteria (we elaborate on the case of a variant of Shannon's entropy) 
therefore yield estimates for the indices. 

We focus here on the case of transient and locally sinusoidal (or harmonic) layers in audio signals, 
using wavelet and local cosine bases. However, the approach we develop may be adapted to different 
signal layers (chirps for example), and to higher dimensions. We provide theoretical estimates for the 
behavior of transientness and tonality indices, and illustrate our results by numerical simulations and tests 
on real sounds. 

II. Theoretical analysis 

We focus on the particular application to audio signals, and limit ourselves to transient and tonal 
features. Our starting point is the assumption that transient signals admit a sparse expansion in a wavelet 
basis (provided the wavelets have small enough time support), and that tonals admit a sparse expansion 
in local cosine basis (with smooth enough window function.) We are naturally led to consider a generic 
redundant "dictionary" made out of two such orthonormal bases, denoted by ip x and w$ respectively (we 
refer to [HI, [HI for detailed tutorials), and signal expansions of the form 

X = ^ OL\ll)\ + J~]p S W S + r = X tr + Xton + T , (1) 

where A and A are (small, and this will be the main sparsity assumption) subsets of the index sets, 
termed significance maps. The nonzero coefficients ct\ are assumed to be independent A/"(0, a x ) random 
variables, and the nonzero coefficients 0s are assumed to be independent A/"(0, as) random variables, r is 
a residual signal, which is not sparse with respect to the two considered bases (we shall talk of spread 
residual), and is to be neglected or described differently. 

Given a signal assumed for simplicity to be of the form (Q3, with unknown values of | A| and |A| (the 
cardinalities of A and A respectively), we are interested in finding estimates for the latter. More precisely, 



we seek estimates for the "transientness" and "tonality" indices 

IM , IA 



hon = taTTTTT • ( 2 ) 



tr |A| + |A| ' ^ ton __ |A| + |A| ' 

We propose a procedure close to the notions of theoretical dimension or ct-entropies, advocated by 

Wickerhauser (see IfTOl for a review.) Our approach is based upon the following heuristics. Consider 

a signal x, and expand it into an orthonormal basis. Estimating the "size" of x in this basis may clearly 

be done by counting the number of nonzero coefficients (the £° norm of the sequence of coefficients), or 

the number of coefficients above some threshold. It has been shown IfTOll that alternative approaches are 

possible, including £ p norms (with p < 2) or entropy, yielding comparable results. Considering a hybrid 

signal as in ([]]) neither its wavelet expansion nor its local cosine expansion will be sparse. However, 

since by assumption only very few coefficients ct\ and (3g are nonzero, most wavelet coefficients (x, i/)\) 

actually originate from the tonal part x ton , and most local cosine coefficients (x, Wg) originate from the 

transient part x tr . Therefore, calculating these £ p norms or entropy from the wavelet coefficients (x, i/)\) is 

expected to provide (approximately) estimates on the number of nonzero (or significant) (3g coefficients, 

and vice versa. We elaborate below on the specific case of the logarithmic dimensions, for which such a 

behavior may be proved. For the sake of simplicity, we shall work in this section in a finite dimensional 

context. 

Definition 1: Given an orthonormal basis B = {e n ,n = 1, . . . N} of a given iV-dimensional signal 

space S = R N , define the logarithmic dimension of x G S in the basis B by 

1 N 
V,3(x) = -Ys^2{\(x,e n )\ 2 ) . (3) 

n=l 

It follows from a simple calculation that in the framework of the signal models under consideration, 

Lemma 1: Given an orthonormal basis B = {e n , n — 1, . . . iV}, assuming that the coefficients (x, e n ) 

are A/"(0, a n ) random variables, one has 

1 - 
V.{V B (x)} = C+-YpM) 1 (4) 

n=l 

where C — 1 +7/ hi(2) & 1.832746177 is a universal constant (7 ~ .5772156649 being Euler's constant.) 

Returning to the model £Q), and assuming that the coefficients a a, A G A and /3g, 5 G A are respectively 

Af(0,a\) and Af(0,ag) independent random variables, the coefficients a\ = (x,i/>\), bs = (x,Wg), are 

zero-mean normal random variables, whose variance depends on whether A G A (or 5 G A) or not. For 



example, 



var{a A } 



^A + E^A^flW'A,^)! 2 if AG A 



E^A^K^A,^)^ 



if Ae'A 



(5) 



and we obtain, for the ^ = {ip\} basis 



E{V*(x)} = C+-\og 2 



n-A+E^i^' 



Ws) 



•AeA 



8eA 



(6) 



\'^A\SeA 

and a similar expression for the logarithmic dimension T>w(x) with respect to the W = {w$} basis. 
In the simpler case where a\ = a, VA G A and a$ = a,V5 E A, we introduce the Parseval weights 



P\ 



(A) = ]T|K^ A >r , p,(A) = £|<ti;,,^ 



(7) 



<5eA 



AeA 



The following property is an immediate consequence of Parseval's formula (i.e. for all /, J^ A | (/, ■?/>> 

•) 



Lemma 2: The Parseval weights satisfy 



< pa(A) < 1 , < p s (A) < 1 . 

Introducing the relative redundancies of the bases \fr and VT with respect to the significance maps 

e(A) = supp A (A) , e(A) = supp 5 (A) , 
AeA <5eA 



(8) 



we obtain simple estimates for the logarithmic dimension. 

Theorem 1: With the above notations, assuming that the coefficients {a\, A G A} and {(3s, 5 G A} 
are independent identically distributed A/"(0, a) and A/"(0, 5") normal variables respectively, and assuming 
r = 0, the following bounds hold 



E{V 9 (x)}>C+^hg 2 (a 2 ) 

+iog 2 (n^V'(A)) 

Vaya 
E{Dv(x)}<C + ^ log 2 (a 2 + e(A)a 2 ) 

+ log 2 (n(^v(A)) 

\A'eA 



1/N 



(9) 



1/7V 



(10) 



with C — 1 + 7/ln(2) ~ 1.832746177. Exchanging the roles of A and A, a similar bound holds for the 

other logarithmic dimension T>w{x). 

Proof: The proposition follows directly from the fact that in such a situation, equation © reduces to 

E{V 9 (x)} = C+\og 2 ( lHa 2 + a 2 p x (A)) 1/N 

v (ID 

xII(-V(A)) 1/iV ), 

from Lemma |2] and the definition of e(A). 6 

This result may be understood and utilized as follows. First notice that the bounds in Equations © 
and (fTOl) differ by |A| log 2 (l + e(A)a 2 /a 2 )/N. Let us assume for a while that this term may be neglected 
(more on that below.) Then the behavior of E {Uy(x)} is essentially controlled by log 2 ( Yl x ,g A (a 2 p\i(A)) 
The behavior of this term is not easy to understand, but a first idea may be obtained by replacing py (A) by 
its "ensemble average" -^ J2\=iP\(A) = jj ^2seA \\ w s\\ 2 = n i which yields the approximate expression: 

E{V 9 (x)}*C + ^log 2 (<T 2 ) 

|A|\ , ^ 2 |A|\ (12) 

- -n) log2 [ a IT) ■ 

Therefore, if the "^/-component" of the signal is sparse enough, i.e. if |A|/iV is sufficiently small (compared 
with 1), E{Vq,(x)} may be expected to behave as log 2 (5" 2 jv)- Set 

Nnf(x) = 2 v ^ x) . (13) 

Replacing V^(x) with its expectation, we see that N^(x) ~ 2 C a 2 ^-, which yields an estimate (up to the 
multiplicative constant 2 c a 2 /N) for the "size" of the tonal component of the signal. Similarly, defining 

N w (x) = 2 Vw{x) (14) 

we obtain a similar estimate (up to the multiplicative constant 2 c a 2 /N) for the "size" of the \I/ component 
of the signal. Both Ny(x) and Nw(%) are computable from the signals wavelet and local cosine expansions, 
and we finally consider their relative proportions, or "rates" 

I - ** f - Jfe Q5) 

J-ton — — , ±tr — — , UJj 

iV* + N w N^ + N w 

which provide the desired estimates for the indices in Eq. ©. 

A few comments are in order here. 



i. The difference between the lower and upper bounds depends on the sparsity |A|/iV of the \&- 
component and the relative redundancy parameters e(A). The latter actually describe the intrinsic 
differences between the two considered bases. When the bases are significantly different, the relative 
redundancy may be expected to be small (notice that in any case, it is smaller than 1.) 
it The relative redundancy parameters e and e differ from the one which is generally considered in 
the literature, namely the coherence M[W U *] = su.p b€Wb , € ^, \(b,b')\ of the dictionary W U ^ 
(see e.g. ifTTTl . [fT2ll . lfT3l .) The latter is intrinsic to the dictionary, while the Parseval weights and 
corresponding e and e provide a finer information, as they also account for the signal models, via 
their dependence in the significance maps A and A. 

Hi. Precise estimates for e and e are difficult to obtain (numerical simulations yield values around 1/4.) 
More precise models for the significance maps A and A could provide better understanding. In 
particular, structured models such as those described in [fi~4l (implementing time persistence in A 
and scale persistence in A) are expected to yield smaller values for the relative redundancies than 
models featuring uniformly distributed significance maps. 
Another interesting point is the sensitivity of such tools with respect to departures from the model, or 

noise. We show that results similar to the above ones still hold true in the presence of white noise, i.e. 

assuming that the residual r in (0Q> is a zero-mean Gaussian white noise. In such a situation, denoting by 

s 2 the variance of the noise r, equation © becomes 



E{Vy(x)} = C+±\og, 



l aga\ <5eA 

x n(E^ 2 K^,^>i 2 +* 2 

A'£A\<5eA 

and a similar expression for the logarithmic dimension T>w{x) with respect to the W = {w$} basis. 
Hence, the approximate expression (fT2l) becomes 

E{V^x)}^C+ l -^\og 2 (a 2 + s 2 ) 

1 |A|\ , / JAI 

The discussion above (suitably adaptated) still holds as long as the signal energy a 2 \ A| exceeds the noise 
energy s 2 N. 

These estimates yield the following algorithm for estimating transientness and tonality indices for sparse 
hybrid signals x as given in ©. 



transientness and tonality indi 




transientness and tonality indices 




Fig. 1. Transientness and tonality estimates for the model (averaged over 10 realizations.) Top: L = 25, and M S {1, . . . 150}; increasing 
curves: I t0 n and Iton', decreasing curves: Itr and Itr- Bottom: M = 25, and L £ {1, . . . 150}; increasing curves: I tr and I tr \ decreasing 
curves: I ton and I to „. 



1) Compute logarithmic dimensions Vq,(x) and V w (x) 

2) Compute N 9 (x) and N w (x) as in £[3]) and (fl4l) . 

3) Compute the estimated rates as in (TTST) . 



III. Numerical results 

We generated several realizations of the signal model (with r = first), with variable numbers L of 
wavelet atoms and fixed number M of local cosines and vice versa, and computed the estimated rates I ton 
and I tr , to be compared with the ground truth ©, i.e. I ton = M/(M + L) and I tr = L/(M + L) = 1 — I t(m . 
As may be seen from Figure \T\ (which corresponds to averages over 10 realizations of the model), the 
estimated curves reproduce quite well the correct ones. Some discrepancies may be observed at the right 
hand side of the curves, where the sparsity assumptions are not valid any more, and the correction terms 
in (fl"2l) come into play. Observe that the curves cross precisely at the correct location M = L. The influence 
of the noise may be seen in Figure |2] a white noise, whose energy equals 30% of the signal's energy, has 
been added. The effect is what can be anticipated from (HO), namely the presence of an additional noise 
term moves the experimental curves away from the theoretical ones. 

Besides the numerical simulations above, the transientness and tonality indices have been tested on real 
audio signals, yielding very sensible results U A first example, based upon a simple ' Castanet' signal (6 s 
long, sampled at 44,100 kHz) is shown in Figure [3] A value for the transientness index and the tonality 
index was computed for all time frames (23 ms long.) Since I ton = 1 — Itr, only the transientness index is 

'Additional material, including sound files, may be found at the web site |http: //www, cmi .univ-mrs ■ f r/~torresan/papers/balance| 



tra nsientness and tonality indices 




tra nsientness and tonality indi 
0.8 




Fig. 2. Influence of white noise: Transientness and tonality estimates for the model (averaged over 10 realizations) with additional white 
noise. Same legends as before. 




Fig. 3. Transientness index for the test 'castanet' signal. Signal (top) and transientness index (bottom.) 



displayed for the sake of clarity. This signal is quite simple, as it essentially exhibits attacks followed by 
harmonic tones, and is thus a "perfect" test for the proposed approach. As may be seen from the bottom 
plot of Figure |3l all attacks are correctly captured, and the corresponding index is quite high. In between 
attacks, the transientness index is very low, which is also natural since the signal is essentially harmonic, 
thus sparsely represented by local cosine basis. 

The second sound example displayed here is a more complex audio signal, extracted from a jazz 
recording (about 6 s. long, sampled at 44,100 kHz) which features "mixed" tonals and transients. The 
numerical results are displayed in Figure @] Notice again that the "obvious" attacks of the signal have been 
captured by the method. A closer examination of the signal (using a "spectrogram type" representation, 
not shown here) shows that in the middle part of the signal (more precisely, between seconds 3 and 5), the 
harmonic content is stronger, which explains the lower average value of I tr there. This illustrates the fact 
that I tr really provides an estimate of the proportion of transients relative to tonals, rather than an absolute 




Fig. 4. Transientness index for the test "jazz" signal. Signal (top) and transientness index (bottom.) 



I 



indicator of the presence of transient, such as the ones used in transient detection [4 J for example! 

More numerical results, in the framework of the hybrid audio coding scheme developed in [3 J, will be 
given in a forthcoming publication |[T4"1 . 



IV. Conclusions 

Sparsity of wavelet and local cosine signal representations may be exploited in order to estimate the 
relative amount of tonal and transient components present in the signal. This approach proves to be 
extremely effective in the context of hybrid audio signal coding 0, [fT5l . and possesses a wider range of 
applications, including image coding fl6]. 

The theoretical analysis we have presented is based on strong a priori assumptions on the signal 
(essentially, a hybrid model such as (Q3, with sparse significance maps A and A.) While this sparsity 
assumption is completely necessary, the equality of variances may be relaxed; in that situation, the indices 
provide estimates on the proportion of energies a 2 |A| and <7 2 |A| of the two layers, rather than their size 
|A| and |A|. This is the case for the numerical results on real signals, for which the variances of the two 
layers are not known. 

Finally, let us simply mention that the approach may be extended to more than two layers, provided 
that the considered orthonormal bases are sufficiently different (in terms of their "Parseval weights", see 
above) to allow the separation. Again, this may prove useful in the context of image coding, where new 
types of waveforms (e.g. curvelets) may be introduced. 

|http : //www, cmi . univ-mrs ■ fr/~t or res an/paper s/balance| . 
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